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ON SOME CLASSES OF BIPARTITE UNITARY OPERATORS 


JULIEN DESCHAMPS, ION NECHITA, AND CLEMENT PELLEGRINI 


Abstract. We investigate unitary operators acting on a tensor product space, with the property 
that the quantum channels they generate, via the Stinespring dilation theorem, are of a particu¬ 
lar type, independently of the state of the ancilla system in the Stinespring relation. The types 
of quantum channels we consider are those of interest in quantum information theory: unitary 
conjugations, constant channels, unital channels, mixed unitary channels, PPT channels, and en¬ 
tanglement breaking channels. For some of the classes of bipartite unitary operators corresponding 
to the above types of channels, we provide explicit characterizations, necessary and/or sufficient 
conditions for membership, and we compute the dimension of the corresponding algebraic variety. 
Inclusions between these classes are considered, and we show that for small dimensions, many of 
these sets are identical. 
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Introduction 

In this work, we study some families of unitary operators acting on a tensor product of two finite 
dimensional Hilbert spaces, having some special properties in relation to the Stinespring dilation 
theorem. This fundamental result in operator algebras [19] states that any linear, completely 
positive, trace preserving map L acting on A4n(C) can be written as 

L(p) = [id®Tr](U(p®/3)[/*), (1) 

where U is a unitary operator acting on the tensor product (S> C^, /3 G A4fc(C) is a positive 
semidefinite matrix of unit trace, and k is a large enough parameter {k = v? suffices). In quantum 
information theory [18], the map L is called a quantum channel, and the matrix /3 is called a density 
matrix (or simply a quantum state). The Hilbert space by which the original space C"' needs 
to be extended is called the environment, or the ancilla space. 

The starting point of our investigation is the remark that the channel L in (1) depends, a priori, 
on the quantum state (5. In the practice of quantum theory, the environment space is usually 
large (most of the times much larger than the system space C”), and thus it is inconvenient to 
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describe the aforementioned dependence of L on /3. More precisely, we would like to characterize 
the unitary operators U, for which, independently on the value of j3, the channel L given by (1) 
belongs to some given class £ of quantum channels. 

In this work, we answer the question above for several classes £ of relevance in quantum in¬ 
formation theory: unitary conjugations V ■ V*, constant channels {L{p) does not depend on the 
quantum state p), unital channels, (£(/) = I), mixed unitary channels (convex combinations of uni¬ 
tary conjugations), PPT channels (channels for which the Choi matrix has a positive semidefinite 
partial transpose), and entanglement breaking channels (channel which, acting on one half of any 
entangled state, yield separable states). In particular within our work, we give a partial positive 
answer to the conjecture formulated in [I, 2]. In [I, 2], it was conjectured that convex combination 
of unitary conjugations can be obtained only by a very specific class of bipartite unitary operators 
(the set Umixed in the present work). We show that this conjecture is true, under some additional 
assumptions. 

It is very important to state at this time that we are not concerned with classes of channels, but 
with classes of bipartite unitary operators. Although these classes are defined in terms of channels, 
we are interested in characterizing the “interaction” unitaries U with the property that, for all 
ancilla states /3, the channel L given by (1) has some fixed set of properties. A similar question was 
studied in [13], where the authors characterize the unitary operators U having the property that 
the only matrices /? which give quantum channels in (1) are quantum states. 

The paper is structured as follows. In Section 1 we define the classes of unitary operators we are 
interested in, and we present some general properties. Sections 2-6 deal each with one or more of 
these classes. In Section 7, we collect some more relations between the different classes, and present 
some equality cases. We close the work with a section containing some open problems. Finally, in 
Appendix A, we discuss the block singular value decomposition of operators. 
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1. Some classes of unitary operators 

Let us fix, once and for all, the space on which the unitary operators we investigate will act. Put 
1-La = Hb = and let Hab = ILa^Hb = €A The compact group of unitary operators 

acting on Hab = will be denoted by Unk or, simply, U: 

U := Ur,k = {C/ G Afnfc(C) I UU* = U*U = I}. 

Starting from the trace linear form Tr : A4fc(C) —)• C and the transposition operation t : 
A4fc(C) —)■ A4fc(C), define the partial trace and the partial transposition respectively by 

TrB = id (g) Tr 
ts = id (g) t. 


In other words 


Tib : Mn{C) ® Mk{C) ^ M„(C) 

{TtB)A 
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and 


ts : Mn{C) ® Mk{C) ^ Mn{C) ® Mk{C) 

A®B^ A®B^. 

Note that we use the following notation for the transposition t{B) = B^. For obvious aesthetic 
reasons, we shall write = ts(X). 

We denote by W1 (C) the set of fc-dimensional density matrices (or quantum states) 

W1^’+(C) := {/? E Affc(C) I /3 > 0 and Tr(/3) = 1}. 

To a unitary transformation U gU and a quantum state /? E Wl|,’^(C), we associate the quantum 
channel : W1„(C) —)• AiniC) defined by the Stinespring formula (see, e.g., [18, Section 8.2.2]) 

LuA^) = Tvb{U-X®/3-U*). ( 2 ) 

Let us introduce some classes of quantum channels having the unitary invariance property: 

VL E £, Fi, ^2 G Un, V2 L{Vi ■ V{)V2* E £. 

We present next a list of such classes of channels, leaving the task of verifying the unitary 
invariance property as an exercise for the reader; 

• Unitary conjugations 

Cant = {X^ VXV*}v^Ur.- 

• Constant channels 

Cconst = {L : Mn{C) -G Mn(C) : 3 a G M^’+(C) s.t. Vp E M^’+(C), L(p) = aj. 

• Unital channels 

Annual = {L : Mn(C) ^ Mn(C) : L{In) = In}. 

• Mixed unitary channels 


I^mixed = convjx hA VXV*]v&j„. 

• Positive partial transpose (PPT) channels 

CppT = {L : Mn{C) ^ Mn{C) : Vp E Mi’+(C), [[L ® id„](p)]^ > 0}. 

• Entanglement breaking channels 

Ceb = {L : MniQ -G Mn(C) : y pG (C), [L (g) id„](p) is separable}. 

We move next to the main definition of this paper, the classes of bipartite unitary channels we 
are interested in. These classes are defined in a natural way as the set of unitary operators inducing, 
via the Stinesping formula (2), independent of the state /3 of the environment system B, quantum 
channels belonging to one of the classes above. This exact notion, in the case of degradable channels 
and entanglement breaking channels, has been considered respectively in [14, Definition 15] and 
[15]. More precisely, we define, for any * E {aut, const, unital, mixed, PPT, EB}, 

U, = {Ug Unk I v/3 E M]:^{C),Lu,p G £4. 


We have, in order: 
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Uaut = {Ue Unk I V/3 G = VfipV^ for some Vp G Un] (3) 

l^const = {U G Unk I V/3 G Lij^p is a constant channel} (4) 

Uumtal = {U ^ Unk \ V/3 G >1 (C) , Lf/,/3(!) = /} (5) 

r 

Urm.ed = {C/ G | V/3 G At(C), Lf/,;,(X) = ^K(/3)C/i(/3)XC/i(/3)* (6) 

withpj(/3) > 0,'Y^pi{/3) = 1 and C/i(/3) £Un} 

i 

UppT = {U e Unk I V/3 G AlV+(C), Lf/,/3 is a PPT channel} (7) 

7/eb = {U G 7/nfc I V/3 G A4],’"''(C), Lf/^^ is an entanglement breaking channel} (8) 

In relation to the class Uaut-, we also define (see Section 2) 

Usingie = {U & Unk \ the Set {Lt/,/3 : /3 G A4^’’^(C)} has 1 element}. (9) 


One of the original motivations of this work was to obtain a characterization of the set Umixed- 
As stepping stones towards a description of this set, we introduce the following classes of bipartite 
unitary operators: 

r 

Uproh = {U^ Unk I 3I/i G Un S.t. V/3 G Mfc’+(C), Lu^X) = Y,PiiP)U^XU* (10) 

2=1 

with pi{(]) > 0 and '^pi{l3) = 1} 

2 

r 

^prob—lin {U G Unk I ^Ui G Un S.t. V/3 G Ml’+iO, Lu^X) = Y^p^U^XU* (11) 

2=1 

with linear functions Pi(/3) > 0 and E PiW) = 1} 

2 

k 

ULck-Mag = {U^Unk\U = YUi® 6*/*, (12) 

2=1 

with Ui G Un and {e*}, {/*} orthonormal bases in C^} 

n 

Ulck-Mag = {U^Unk\U = Y ^ift ® U,, (13) 

2=1 

with Ui G Uk and {cj}, {/*} orthonormal bases in C”} 

Let us mention at this point a simple but fundamental property of the sets of bipartite unitary 
matrices we have just introduced. 

Lemma 1.1. The local unitary group Un x Uk acts by left and right multiplication on U^, for all 
* G {aut, const, unital, mixed, PPT, EB, single, prob,prob — lin, block — diag — A, block — diag — B}: 

yu €U„'dVi,V 2 £Un,yWi,W 2 £Uk, iVi 0 Wi)U{V 2 ®W 2 ) £U,. 

Proof. In the case of * G {aut, const, unital, mixed, PPT, EB}, the proof follows from the bi¬ 
unitary invariance of the corresponding class E* and from the fact that in the definition of 7/*, we 
require the condition to hold for all states on the environment /3 G Al|.’'*~(C). The other cases are 
easy verifications, we leave the details to the reader. □ 
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As suggested by the property above, the subgroup Un®hlk ^ Unk plays an important role in our 
study. It turns out that the flip operator 

: C*" ® C*" ^ C" 0 C^, FnX®y = y®x (14) 

will also be of particular importance, in light of the following classical theorem. 

Theorem 1.2. (^[8, Theorem 3.1]J Let G be a compact group such thatlAn ®tlk^G F lA^k- Then 
G is one of the following 

(1) G = lAn 0 Gk; 

(2) G = Unki 

(3) Ifn = k,G={Un®Un,Fn). 


By direct computation, one can show that the following chain of inclusions holds : 


U, 


C C u 


block—diag — ^prob — ^mixed — ^unital — 




CU. 


Note that one can define “i?”-versions of the above sets, in an obvious way, by swapping the 
tensor factors A and B ( above inclusions are still true for U^). 

One of the main focuses of the current work will be to understand which are the inclusions above 
which are strict and which are actually equalities. 

We end this section by showing that if two interaction unitary operators U, V generate the same 
channels for all states /? on the environment, then they are related by a unitary operator acting on 
the environment C^; this result will turn out to be useful later on. 


Lemma 1.3. Let U,V € Unk be two bipartite unitary operators. The following two assertions are 
equivalent: 

(1) For all density matrix on ft ^ the quantum channels U and V induce are equal: 

LuAp) = LvAp) /o’" P ^ A4n’’^(C); 

(2) There exists a unitary operator W £Uk such that U = {I ® W)V. 

Proof. We only prove “1 2”, since the converse follows from direct calculation. We start form 

the hypothesis, 

Vp G Mn’+(C), G [id 0 TV](C/p 0/31/*) = [id 0 Tr](yp 0/3I/*). 

By linearity, we can replace p, resp. ft, by arbitrary complex matrices A G A4n(C), resp. B G 
A4fc(C). Again using linearity, we replace the simple tensor A ® B by a general element X G 
M.n{C) 0 Alfc(C) to obtain the first line below, and we continue by obvious equivalent reformula¬ 
tions: 


VXGMnk(C), 
VX G Mnk(C), VA G Mn(C), 
yx G Mnk(C), VA G Mn(C), 
VA G Mn(C), 
VA G Mn(C), 
3WeMk(C), 


[id0Tr](C/AC/*) = [id 0 TV](LAL*) 
Tr(C/XC/*A 0 I) = Tt:{VXV*A 0 I) 
Tr{XU*A 0 lU) = Tr(AWA 0 IV) 
U*A®IU = V*A®IV 
A®I{UV*) = {UV*)A®I 
uv* = I ®w 


and we conclude the proof by noticing that W has to be unitary, since UV* is. 


□ 


2. Bipartite unitary operators producing unitary conjugations and constant 

CHANNELS 

In this section we provide characterizations of the sets Uaut (3), Ugingie (9); and lAconst (4), 
showing that only tensor products of unitary operators (resp. flipped tensor products) belong to 
these classes. 
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We start with the set of automorphisms, that is the set of unitary conjugation channels. 


Theorem 2.1. Let U G U^k « bipartite unitary operator such that, for all quantum states /3 G 
there exists an unitary operator Vjs £ Lin such that Lu^p{p) = VppV^ for all p G A4n~^(C). 
Then, there exist unitary operators V £Un and W £Uk such that U = V . In other words, 

Uaut = {V®W : V £ Un, G Uk}. 


Proof. The proof is easy, and consists of two steps: we show first that the unitary operators Vp can 
be chosen to not depend on /3, and then, using Lemma 1.3, we show that the unitary U has the 
required form. 

Let us first introduce some notation. To any matrix 

n 

Mn(C) 3 A = Aijae*, 
ij=l 

associate its vectorization a = vect(^) defined by 

n 

C"" (8> C"" 9 a = AijCi 0 ej, 
hi=i 

where is some fixed orthonormal basis of C”. 

Denote V := 1//^, the unitary which appears for (3 = I/k, and let v = vect(D). For an arbitrary 
(3 £ let /3 = {/3 +1/k)/2. Since quantum channels are linear in (3, we have that 

Vp G Mi’+(C), 2VppV^ = V-^pVi + VpV*. 

At the level of Choi matrices, the above equation reads 


2vpv*p = v-pv*^ + vv*. 

On the left side of the equation above the operator is a rank one operator. This way, in order that 
the sum in the right side is a rank one operator, the vector has to be proportional. Since all the 
involved vectors are of norm ^/n, it follows that they should all be the same, up to a phase. The 
same holds for the unitary operators, which concludes the first step of the proof. 

Let U := V 0 I. It is easy to see that 

Vp G Mi’+(C), V/3 G M^’+(C), L^^^(p) = VpV* = LuAp), 


so, by Lemma 1.3, there exists an unitary operator W £ Uk such that U = {I ® W)U, and the 
proof is complete. □ 


We show next that the class Ugingie, he. the set of unitary operators U £ Unk with the property 
that the map (3 i-a Ljj^p is constant, is actually identical to Uaut- 
For the proof of Theorem 2.3, we need the following lemma. 


Lemma 2.2. Let M £ A4n(C)(8> Alfc(C) be a given matrix. The following conditions are equivalent: 

(1) For all matrices X £ M.n{C) and Y £ Alfc(C) such that TrT = 0, we have Tic{M ■ X <SiY) = 

0; 

(2) There exist a matrix A £ M.n{C) such that M = A<Si Ik- 
Proof. The non-trivial implication follows from the following equation 

[Mn{C) ® (Mfc(C) 0 C4)]^ = Mn{C) 0 Clk. 

□ 
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Theorem 2.3. Let U G Unk o bipartite unitary operator such that, for all quantum states p G 
7Wn''~(C) and /?, 7 G we have that Lu^pip) = Lu^^{p). Then, there exist unitary operators 

V G ILn and W GUk such that U = V . In other words, the set ILsingie defined in (9) is equal 
to 

Usingle = {V ®W \ V G hin, W G Uk} = Kaut- 
Proof. By using linearity, the hypothesis translates to the following equality 

VX, Y G Mn{C), VZ G Mk{C) s.t. TvZ = 0, Tv[U{X ® Z)U*{Y ® 4)] = 0. 

By reshaping the operator U, the previous equality can be re-written as 

Tv[Y^ ^X^Z-11*11] =0. (15) 

where U G Ai^xn^ki^) is the operator 

n k 

u=Y. E ® ® ® ® 4*)’ 

*7=1 a;,y=l 

see Figure 1 for a graphical representation of U and of the previous equality (for the Penrose 
graphical notation for tensors, in the form used here, see [ 6 , Section 3]). From Lemma 2.2, it 
follows that there exist an operator A G Af„ 2 (C) such that U*U = A® Ik- Since the rank of U*U 
is at most k, it follows that A has rank at most 1, i.e. there exist a vector a G C” ® C” such that 
A = aa*. From the equality U*U = aa* ® Ik, we deduce that there exists an operator W GUk such 
that JJ = a* ®W . Writing 

n 

Mn{C)3V= {Cj ® ej,a)eie*, 
i,j=^ 

we get U = V ® W. The fact that V GUn follows from the unitarity of U , and this concludes the 
proof of the first implication. The fact that tensor product unitary operators belong to Ugingie 
be verified by direct computation. □ 



Figure 1. The reshaping [7 of a unitary operator U (left) and the trace inequality (right). 


For the case of Uconst, it is easy to see that it is related to the previous case via a flip operation, 
although there is a slight technical complication. This question has also been considered in [15]. 

Theorem 2.4. Consider the set Uconst of bipartite unitary operators such that, for all quantum 
states p,cr G {C) and fi G we have that Ljj^p^p) = Ljj^p{a). 

If k rn for r = 1,2,..., then Uconst is empty. 

If k = rn for some positive r, then we have 

Uconst — {{In ® V){Fn ® Ir){In ® IF) : F, W G Uk), (16) 

where Fn G Un^ denotes the flip operator (14); see Figure 2 for the diagrammatic representation of 
such an operator. 
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• D 

u 

♦ » 


Figure 2. The diagrammatic representation of a general element from Uconst, in 
the case where k = nr. The diamond shaped labels correspond to the vector space 
CT 



Proof. Let us start with the easy implication, considering an operator U as in (16). By direct 
computation, one can see that 

Lu,p{p) = [idn ® Tr,.](/3), 

proving that, for arbitrary /3, the channel is constant. 

The proof of the difficult implication starts in the same way as the one of Theorem 2.3. The 
hypothesis that the channels Lu^js are constant translates to the following condition: 

VX,y e Mn(C), VZ e Mk{C) s.t. TrX = 0, Tr[U{X^Z)U*{Y(^Ik)] = 0. 


As before, after reshaping the matrix U into U, the previous relation becomes (15), with the 
difference that this time, we have TrA = 0. Using Lemma 2.2, we conclude that there exists an 
operator A G Af„fe(C) such that 

[>*[> = 42 ) (17) 

where the superscripts indicate on which factor of the tensor products the operators are acting 
{A acts on the first copy of C"’ and on while the identity operator acts on the second copy of 
C"'). Using the fact that U is unitary, we have that the rank of the matrix U*U is precisely k, so 
we must have k = n x rank(A). We conclude that if k is not a multiple of n, then operators U 
with the above property cannot exist. We put now r := rank(A), so that k = rn. Since U*U is a 
positive semidefinite matrix, A is a positive semidefinite matrix of rank r, so we write A = B*B, 
for a matrix B G Xirxnk(C) = Ai^xn^ri'^)- The fact that U is unitary translates to the following 
equality 

BB* = Inr, (18) 

where 

n r 

Mnr{C) 3 ^ ^ {g^, Bci ® Cj ® gy)ei ® gx ■ e* ® 5 *, 

ij=l x,y=l 

for some orthonormal basis {gi,... gr} of C'"; see Figure 3 for a graphical representation of the 
previous equalities. Thus W := B is a unitary operator. From equation (17), we find that there is 
another unitary operaor V G Unr such that 


[/(h2,3) 



where the three tensor legs correspond to C” ® = C”" ® C"" ® C’’. 


The conclusion follows. 


□ 



= A 


Figure 3. A diagramatic representation of (18), from which the unitarity of B follows. 
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Remark 2.5. If n = k in the previous result, equation (16) can he written as U = Fn{V 0 W) for 
a pair of unitary operators V,W Glln, so in this case we have 

Hconst — FjiLIsingle — hlgingleFn- 


3. Bipartite unitary operators producing unital channels 

In this section we study the set Uunitai bipartite unitary operations which yield unital channels 
for every choice of the state on the auxiliary space. 

Using linearity, one can extend the definition (5) to the whole space of A: x A: complex matrices: 

Uunitai = {UGU\yBG Mk(C), Lu,b(I) = (TrB)!}. 


Theorem 3.1. One has 

In paHicular, 


h^unital hifik ^nk' 


Proof. Let us first show the “C” inclusion in the above equality. Take U G Hunitai and put V = 
G and W = VV*. One has, for any B G Alfc(C), 


{Ti-B)In = Lu,B{In) = [id ® Tr](U • I ® B • U*) 
= [id®Tr](UU* 

= [id®Tr](IU-I®R^) 


By block-decomposing W in an arbitrary orthonormal basis e* of 

k 


i,j=i 


we get that for all B G Alfe(C), 


k 

{TtB)In = ^ Wij ■ {ej,Bei). 
hi=i 


(19) 

( 20 ) 


Choosing B = CiC*, we get Wij = Sijin and hence W = Ink- In other words, V = G link, which 
finishes the proof of the first inclusion. 

The second inclusion follows by working backwards the previous arguments: since V = G link, 
equations (20) and (19) hold. □ 


Since both sets link and ll^p. are algebraic varieties (i.e. they can be describes as the zero-set of 
a system of polynomial equations), we obtain the following corollary. 


Corollary 3.2. The setUunitai is a real algebraic variety. 


We are going to investigate next Uunitai, as an algebraic variety. We compute first the dimension 
of the “enveloping tangent space” of Uunitai at a point which is a block-diagonal unitary 


k 

U = Y.Ui®e,f* 

i=\ 


The notion of enveloping tangent spaces was introduced in [20] (also called defect), and it is simply 
defined by (see also [4]) 


Tu{Uunitai) ■= Tu{Unk) C TuiU^/.). 
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Proposition 3.3. The dimension of the enveloping tangent space ofUunitai cit a point which is a 
block-diagonal unitary of the form 

k 

U = Y,Ui®eif* 

i=l 


is given by 


Du = ^ dl, 


^i3 1 


( 21 ) 


where Aij is the set {{Xx,dx)} where \x is an eigenvalue of the unitary operator UiU* having 
multiplicity dx- 


Proof. A matrix A G J\Ank(C) is an element of the enveloping tangent space at U if and only if 
both matrices U + eA and (U + eA)^ are unitary, up to the first order in e. The unitarity of U -j-sA 
is equivalent to the condition UA* + AU* = 0, while the unitarity of {U + eA)^ is equivalent to 
U{A^)* + A^U* = 0 (note that we have used the fact U = U^). 

Writing A as a block matrix 

^ ® e*/*, 

ij 

the first condition C/A* + AU* = 0 is equivalent to the following system of equations 

Vi,j, UiA*, + Ai,U* = 0, 
while the condition C/(A^)* + Al'U* = 0 is equivalent to 

VC,j, UiA*^ + A,iU* = U 

First, let us note that the diagonal blocks An appear only in two identical equations 

UA*, + AnU*=U 

The general solution to the equation above is An = BnUi, where Bn is an arbitrary anti-hermitian 
matrix (Bn + B*- = 0). Hence, the total dimension of the diagonal blocks of A is kn^. Note that 
this corresponds to the case i = j \n formula (21): there, An = {(l,n)}. 

Let us now study off-diagonal blocks of A. Again, the equations are decoupled: for i < j, one 
has to solve 


UiA*, + AijU*=t) ( 22 ) 

C/iAL + Aj^U* = 0. (23) 

From the first equation, one finds Aji = —UjA*-Ui. Plugging this into the second equation, we 
have to solve now 

RAj^ j — AijS = 0, 

where R = C/jC// and S = U*Ui. This is the well-known Sylvester equation. From the analysis in 
[10, Chapter VIII], the dimension of the solution space of this homogenous equation depends of the 
Jordan block structure of the matrices R and S. Since in our case both R and S are unitary (hence 
diagonalizable), the Jordan blocks have unit dimension. Moreover, R and S have the same spectrum 
Aij. It follows from [10, Chapter VIII, eq. (19)] that the complex dimension of the solutions of the 
system (22)-(23) is precisely 

( Aaj ^dx )6Ay 

and the proof is complete. □ 

The proof above can be adapted mutatis mutandis to the case of (i?)-classical unitary operators, 
as follows. 
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Corollary 3.4. The same result holds for B-block-diagonal unitary operators of the from 

n 

U = Y,eif*®Ui. 

i=l 

Corollary 3.5. The dimension of the enveloping tangent space of Uunitai o product unitary 
operator U = V ® W is , which is also the dimension ofUnk- 

Corollary 3.6. For k = 2, the dimension of the enveloping tangent space of Uunitai at a point 
U = I (BV = I ®eif^ + V ® 62 / 2 * is 

Di®v = 2 n^ + 2 ^ d\, 

A 

where d\ are the multiplicities of the eigenvalues A ofV. 

Corollary 3.7. Consider a block-diagonal unitary operator 

k 

U = Y,Ui®eift, 

i=l 

where the operators Ui are in generic position; 

Vi / j, has a simple spectrum. 

The dimension of the enveloping tangent space of Uunitai at U is then 

Du = kn^ + — nk. (24) 

Note that the expression above is symmetric in n and k. 

We conjecture that the expression (24) is the dimension of UunitaU a-s an algebraic variety, see 
Conjecture 8.1. 


4. Bipartite unitary operators producing PPT channels 

We consider in this section PPT channels and bipartite unitary operators which produce such 
channels via the Stinespring formula, independent of the state of the environment. 

Recall that the maximally entangled state is the matrix (here, we drop the normalization con¬ 
stant) 

n 

A4„2(C) 3 := ^ aej* ® Cicj*. 

A quantum channel L is said to be PPT if and only if its Choi matrix 

Cl ■■= [id® L](nn) 

is PPT, i.e. C^ > 0. Hence, the set Uppt admits the following characterization: 

UppT = {U G Unk ■ the map (3 1 —>■ C^^^ is positive}. 



Figure 4. From left to right: the matrix the Choi matrix of the map /3 1 — )• 

and the same Choi matrix, with U replaced by . 
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Since the structure of positive maps between matrix algebras is rather poorly understood, we 
focus for the moment on a subset oiUppT, namely 

UppT ^ ^CPPT ■= {U G Unk ■ the map (3 i—)• Cp^ ^ is completely positive}. 

We have the following description of the set Ucppt, in which, remarkably, the partial transpose 
of U plays a special role. 

Proposition 4.1. For all n,k, we have 

Ucppt = {U G Unk : (4 ® U^){Fn ® 4)(4 ® U^T > Oj- 

Proof. We use again the fact that complete positivity is characterized by the fact that the Choi 
matrix is positive semidefinite. In Figure 4, we have depicted in the left image the matrix 
while in the center panel we have the Choi matrix of the map ft i—>■ Cp^^. The right-most panel 
contains the diagram of the same Choi matrix, where we have replaced U by its partial transpose 
, in order to obtain a nicer expression. The equality of the last two panels contains the proof of 
the claim. □ 

In order to further simplify the description given above, by conjugating the above expression by 
the pseudo-inverse of the matrix U^, we are focusing next on the study of the set 

VcppT := {P G Mnt(C) orthogonal projection : (4 ® P)(4i ® 4)(4 0 P) > 0 }, (25) 

and we have Ucppt = {U G Unk '■ Ran([/^) G Pcppt}- 

We have gathered the following properties of the set Pcppt] we leave the proofs of these simple 
facts to the reader. 

(1) It is locally unitarily invariant: for all U GUn and V G Uk, if P = {U 0 V)P{U 0 V)* , then 
(4 ® P)iFn 0 4)(4 0P) = {U0U0V){In0 P){Fn 0 hWn 0P){U 0U 0V)*. 

(2) Ink ^ Pcppt- In other words, no product unitary lies inside Ucppt, nor inside Uppp- As 
a consequence, we have Uunitai fiUppp = 0 - 

(3) Since Uamst P Uppp, if k = nr, for any unitary operator V G Unr, Pv £ Pcppt, where Py 
is depicted in Figure 5. 


• 

- . 

_ 

1 • 


4 

0 

n 

V 


V* 

♦ 


> - ♦ 

:> 


Figure 5. For any unitary operator V G Unr, the orthogonal projection Py depicted 
here is an element of Pcppt- 

(4) \f x0 y G RanP, then, for any x' G C" such that x' F x, x' 0 y ^ RanP. 

(5) For any x G C", ||a:|| = I, and any orthogonal projection Q G M.k{C), xx* 0 Q G PcPPT- 
At the level of examples, the only observation here is that Uconst ff Uppp- We refer the reader 

to Section 8 for some related open problems. 

5. Bipartite unitary operators producing mixed unitary channels 

In this section we investigate the set Umixed- We provide necessary conditions for a bipartite 
unitary operator U to belong to Umixed, and we show that in the case of qubit channels (n = 2), 
the sets Umixed and Uunitai are equal. 

Recall that the Kraus operator space of a quantum channel L{p) = Yli PiPPt is the space 
K{L) = span{4} [9, 22]; note that K[L) does not depend on the choice of Kraus operators for L, 
since all Kraus representations are related by unitary transformations [18, Theorem 8.2]. One of 
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the main observations in [22] was that for a mixed unitary channel L, K{L) nUn ^ The next 
result builds on this remark. 

Proposition 5.1. Let U G Umixed be a bipartite unitary operator. Then, for any unit vector f G 
and any orthonormal basis {eff of , we have 

span {In ® e* ■ U ■ In ^ /}^=i n / 0- 

Proof. For any choice of / and {cj}, the operators 

Ei\=In®e*-U-In®f 

are Kraus operators for the channel Ljjjf*. Since the channel is mixed unitary, it follows from [22, 
Section IV] that the linear span of the Ei should contain a unitary operator. □ 

Remark 5.2. Note that in the statement above, the set 

span {In ®e* - U ■ In® f}i^i 

does not depend on the particular choice of the basis {ci}, but only on the vector f. 

As a direct consequence of the above result, we obtain the following simple criterion for deciding 
if a given unitary matrix U is an element of Umixed- 


Corollary 5.3. Let U G lAnk be a bipartite unitary operator with the following property: 

3/ G s.h Ve G C^ In® e* ■ U ■ In® f ^Un- 

Then, U Umixed' 


With the help of the criterion above, we present next an example of an element U G Ll^iock-diag \ 
Umixedi which shows, in particular, that the inclusion Umixed ® Uunitai is strict; this example is 
motivated by [16, Section 4.3] and [22, Example 1]. Let U G 11^.2 be 


0 l' 

, 1 

'l 1 ■ 

* 1 

1 i 

1 0 

+ ^3^3®^ 

1 -1 

+ 6464 (g) —= 

V2 

-i -1 


U = eie^ ®I2 + 626*2 

Obviously, U G U^. In the spirit of the criterion above, compute 

/ z + w z — iw 
= Diag z, w, 


U® \z w'\ - U ■ li' 




V2 ’ v/2 


Asking for the diagonal matrix above to be unitary leads to a contradiction, and thus, by Corollary 
5.3, we conclude U ^ Umixed- 


Let us now consider the qubit case n = 2, which is special because the quantum Birkhoff result 
holds for qubits [21, 16, 17]. 

Proposition 5.4. Let L : A42(C) —)• M. 2 {C) be a completely positive, unital and trace preserving 
map. Then, there exist unitary operators Ui,... ,Uk £Un and probabilities pi > 0, YliPi — ^ such 
that 

k 

VXeM2(C), L{X) = Y,PiUiXU*. 

i=l 

As a corollary, since every unital channel Ljj^jj must be mixed unitary, we obtain the following 
result. 


Proposition 5.5. For n = 2 and any k >2, we have that 

Umixed Unnital- 

In particular. U G Umixed iff U G ZY n U^. 
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6. Block-diagonal bipartite unitary operations 

In this section we study the set of block-diagonal operators, ^biock-diag- Before proving any 
results on this class, let us provide another way of writing equation ( 12 ), which has the benefit of 
being unique in a certain sense. As a corollary, we deduce that the only unitary transformations 
which are blcok-diagonal with respect to both sub-systems A and B are given by partial isometries. 


Definition 6.1. Two unitary operators U and V are said to be in relation, denoted by U y, if 
there exists a eonstant \ in C with |A| = 1, sueh that 

U = XV. 


Proposition 6.2. A bipartite unitary transformation U € U is an element of if and 

only if it ean be written as 

r 

U = Y,Ui®Ru (26) 

i=l 

where Ui are unitary operators aeting on C” and Ri are partial isometries Ri : ^ £f sueh that 

r r 

Y,R,R* = Y,R*Ri = h (27) 

i=l i=l 

and Ui oo Uj for all i ^ j. 

Moreover, the decomposition (26) is unique, up to ~ and permutation of the terms in the sum. 
Proof. Consider two decompositions of a same operator in form of (26) 

r s 

'^Ui® Ri = '^Vj® Qj, (28) 

i=l j=l 


with 

r r s s 

Y, RiRi = Y = E = E ’ (29) 

Z =1 2=1 j = l j = l 

Up oo Uq for all p 7 ^ g and V) I 4 n for all I m. 

For all z in {1,..., r} and j in {1 ,. .., s}, applying I ® R* on the left and I ® Q*jQj on the right. 
Equation (28) becomes 


Ui ® R*iRiQ*Qj = Vj ® R*QjQ*Qj = Vj ® R*Qj . 
This particularly means that 


R*Qj = 0 


or 


3Xij, |Ajj| — 1, 


Ui — Xij Vj 


and R*RiQ*Qj = ^R*Q 


X 




Now note that we have 

S 

Ui ® R*Ri = YUi® R*iRiQ*Qj + 0 . 

This implies that at least one of the terms in the sum is non-trivial. Moreover, since V) Vm for 
all I 7 ^ m, the operator Ui can be in relation with only one of the Vj's. Therefore, we obtain r = s 
and for all i, there exist a unique j such that Ui = XijVj and R*RiQ*Qj = 1/Xij R*Qj. After 
following the same strategy with I ® QjQ* on the left and I ® R* on the right, we now can deduce 
that Ri = IjXijQj. The result follows □ 


Another point of view on block-diagonal unitaries is the fact captured in the next proposition. 
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Proposition 6.3. A bipartite unitary operation is block-diagonal if and only if it admits a block- 
singular value decomposition with respect to B: 

^block—diag ^nk block—diag^ 

where is the set (see Appendix A) 

k 

•^^block—diag {^ = E Xi <8) Cif*, for Xi G Atn(C) and orthonormal bases ei, fi of C^}. 
i=l 

In particular, the sets ^Sock-diag> ^ are algebraic varieties. 

Proof For U G M^i^^^-diag^ write 

k 

UU* = Y,XiX* ^ae*. (31) 

i=l 

The above matrix is the identity if and only if each of its diagonal blocks XiX* is the identity, and 
the claim follows. □ 


Let us now investigate the relation between the two classes ^^ocfc-dmg ^Uock-diag- start 
by presenting an algorithm allowing to check if a unitary matrix U in Unk belongs to hiuock-diag- 
This key result relies on Theorem A.l. 


Proposition 6.4. Let U be in Unk- Consider the operators {Xa)^^i n^ in Alfc(C) defined, for 
all a = 1, , by 

;= [TV®id]((S* ®/) U), 

with an orthonormal basis of M.n{C). Then, U belongs to ^^ocfc-dmg if o.'^d only if the 

2 2 

families resp. {Xf^Xp}f^ consist of commuting, normal operators. 

Proof. Fix t/ G Unk- Thanks to Theorem A.l, the matrix U has a block-diagonal SVD, that is, there 
exists matrices Ui, .. .Up G M.n{C) and partial isometries Ri, ... ,Rp G Alfe(C) having orthogonal 
initial, resp. final, projections such that 

p 

U = ^Ui®Ri, (32) 

i=l 

2 2 

if and only if the families resp. consist of commuting, normal opera¬ 
tors. Then the unitarity of the (Lj)i=i,,,.,p’s directly follows from the unitarity of U. □ 

As a first application of the above result let us present the relation between 
^biock-diag ™ qubit space C^. 

Proposition 6.5. If n = 2, then 

(- ijA 

^block—diag — ^block—diag' 

Proof. Any element U G written as 

U = eifl ®Ui + 62/2 ® U2. 

In order to apply Proposition 6.4 we consider the orthonormal basis {Eij = fie*,i,j = 1,2}. In 
particular 

Xij = dijUj, 
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for all i,j = 1,2. Consider now the sets 

i,j, k,l = 1,2} = {I, UiU^,U 2 Un and 

{X:^Xki,i,j,k,l = 1,2} = {I,U}U2,U^Ui}. 

It is obvious that these sets consist of commuting, normal operators, finishing the proof. □ 

Corollary 6.6. In the case n = k = 2, we have 

^block—diag ^block—diag' 

Note however that the inclusion in the above result is strict (in the case n > 3, k = 2). For 
n = 3, and arbitrary k > 2, we construct next an example of a unitary operator being in hi^ock-diag 
but not in U^iock-diag- 

Consider two non commuting unitary operators V, W in A^fc(C) and an orthonormal basis (cj) 
of and define 

[/ = eie} 0 / + eae^ 0 C + 636 ^ 0 W. 

By construction, the operator U belongs to check that it is not in ^^ocfc-diag- 

Consider the orthonormal basis {Eij = eie},i,j = 1, 2, 3} of Ai^^C), we have for example 

Xii = I, X 22 = V, X 33 = VF. 

We immediately note that the operators Xl]^X 22 = V and X};^X 33 = VF do not commute. Since 
this commutativity is necessary to be in U^ock-diag (Proposition 6.4), we conclude that U doesn’t 
belong to 

Another class of interesting block-diagonal (with respect to the second system, B) operators are 
circulant unitary matrices. 

Proposition 6.7. Let X G M.nk{E.) he a circulant matrix. Then X G -^hiock-diag- particular, 
any circulant unitary operator U is block-diagonal with respect to the B factor. 

Proof. In the proof of this result, since we are going to make us of circularity properties, the indices 
for the matrices we consider are starting at zero. Recall that a matrix X G Ainki^-) is circulant iff 

V 0 < i, j < nk, Xij = , 

where x G is the first row of X and we write [a]p = a mod p, for any integers a and p. 
Circulant matrices are known to be precisely the matrices which are diagonal in the Fourier basis. 
Recall that the Fourier matrix (which implements the change of bases between the canonical basis 
and the Fourier basis) is given by 

yO<i,j<p, Fp{i,j)=uj^\ 

where uj = exp(27rz/p) is a primitive p-th root of unity. Finally, since we are going to work 
with matrices living in a tensor product space, the element (s,t) of the block (f,j) of a matrix 
A G A4ji(C) 0 A4fc(C) is Ai.k+s,j-k-{-t- 

Now that the notation is fixed, consider a circular unitary matrix X G Aink(C) and let x be its 
first row. For any matrix A G Aik{C), define 

Aa = [id0T4](A-/0 A). 

We show next that the matrices X^ are all circulant, fact which, by Theorem A.l, suffices to 
conclude, since all the matrices appearing in the theorem will be simultaneously diagonalizable in 
the Fourier basis. 
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For all 0 < z, j < n, we have 


fc-i 

s,t=0 


k-1 

s,t=0 


The crucial observation is that the above quantity only depends on the difference j — i- indeed, if 
[{j — i) = {j' — then there exists some r such that j' — i' = j — i + nr, and thus 

[(/ - i')k + {t- s)]nk = [{j - i)k + {t- s)+ nkr]nk = [if - i')k + {t - s)]nk, 


showing that the matrix Xj^ is circular, and finishing the proof. 

The statement about circular unitary operators follows from the general case using Proposition 


6.3. 


□ 


We turn next to the study of the unitary operators which are block-diagonal with respect to 
both systems A and B. 

Proposition 6.8. A unitary operator U is block diagonal with respect to both tensor factors A and 
B ft.e. U e iff 

s r 

U = ^ ^ ^ ^ ^ij Qi ® Rj , 

i=l j=l 

where, for all i = 1,... ,s, j = 1,... ,r, |Ajj| = 1, and where iQi)i=i^...^s , iRj)j=i,...,r are two family 
of orthogonal partial isometries respectively on CA and satisfying 

s s r r 

= = = (33) 

1=1 i=l j=l j=l 

Proof. Let U be an element in the intersection ht^ock-diag ^ kl^iock-diag- Then, U admits both 
decompositions 

s r 

U = '^Qi®Vi = '^Uj®Rj. 
i=l j=l 

Applying Q* ® Rt on the left, we obtain 

Q*Qi 0 R*Vi = Q*Uj 0 R*Rj . 

Then there exists fiij 0 such that 

Q*iQi = IPjQ*i Uj and R*Vi = — R*Rj . 

Hij 

Now since the Qi'^s and the Rj^s satisfy (33) we end up with 

s 1 

Uj — ^ ^ f^ij Qi and ^ ^ Rj • 

1=1 j=i 

Since the operators Uj and V) are unitary, we conclude that = 1 and that gives the result. □ 
Finally, we compute next the (real) dimension and Uj^ck-diag^Kiock-diag- 
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Proposition 6.9. The real dimension of the algebraic variety 

..A ifn = l 

^^^Uuock-diag I ^^^2 ^ 2k - 2) if n > 1. 

The real dimension of the algebraic variety 

= 1 
> 1 . 


{ {nk)"^ 

2n2 + 2k‘^ + nk-2n-2k 


if min(n, k) 
if min(n, k) 


Proof. Let us first perform a heuristic parameter counting for a generic element 

k 

ULck-diag^U = Y,Ui®eifl 
i=l 

The choice of the two orthonormal bases {e*} and {fi} in (32) corresponds to a total of 2k‘^ real 
parameters, since dimji^fc = Each matrix Ui accounts for real parameters, so, in total, we 
get extra real parameters. However, in each term Ui (8> eiff, two of the three complex phases 
of Xi,ei, fi are redundant, so we have over counted 2k real parameters. We conclude that the real 
dimension of should be 2k‘^ + kn^ — 2k. The above reasoning over-counts in the case 

n = 1. Indeed, in that case one can ignore the matrices Ui (the phase can be absorbed in the Ciff 
part), and we are left with an unitary matrix, which counts for k"^ real parameters. The rigourous 
proof of this result is very similar to the one of Proposition A. 5, and is left to the reader. 

Let us now find the dimension of the intersection. Similarly, let us count parameters for a generic 
element of the form 

n k 

^Uock-diag ^block-diag ^ ^ ^ X/ X/ ^T®*/** ® Sjh*j 

i=l j=l 

where, for all i = l,...,n, j = 1,...,A;, |Ajj| = 1, and (e*), {fi), {gj) and (hj) are orthonormal 
bases of C"’ and C^, respectively. 

The choice of the four orthonormal bases corresponds to a total of 2n‘^ + 2k‘^ real parameters, the 
choice of the coefficients nk = 2nk — nk. Since, in XijCiff 0 gjht, all the phases can be absorbed 
in the coefficient \ij, we have over counted 2n + 2k real parameters. Again, the case min(n, A;) = 1 
is degenerated, since any unitary operator is of the desired form. □ 


7. Further relations between unitary classes 


As discussed in the introduction, the following chain of inclusions holds: 

^block—diag — ^probl—lin U hiprob ® himixed U Idunital- (^‘^) 

We discuss in this section the situations when some of the above inclusions are equalities. See 
Section 8 for some related open questions. 


Theorem 7.1. The sets U^ock-diag ^prob-iin are equal. 

Proof. Without loss of generality, we assume that the unitary operators Ui are different up to a 
phase, i.e. Ui ^ Uj, for all i / j. 

Due to the decomposition of any matrix into a linear combination of at most four matrices in 
Al|,’''~(C) (positive and negative hermitian parts and their equivalents for the anti-hermitian part), 
the Pi’s can be extended by linearity to positive functionals on Alfc(C). By Riesz Theorem, for 
each i, there exists a matrix Mi such that 


VA G Mfc(C), p^{X) = TriMiX). 
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Note now that the values Pi(/3) = Tr(Mj/?) are non-negative for all /3 in Therefore the 

matrix Mj is actually positive semi-definite since we have 

Tr(Mi/3) = pi{P) = p^if) = Tr(Mi/3) = Tr(M;/3), V/3 G >1^’+(C). 

We then conclude that the matrices Mi can be written as Mi = RiR* = R*iRi = Rf for some 
hermitian and positive semi-definite Ri. 

Moreover, using that for all j3 in 


i=l 


[E«* 


2 = 1 


= 1 , 


it follows that ~ Then, for all (3 in we can write the quantum channel Lu^j^ 

as 

r 

LuAX) = TVb(C/ • X ® /3 • [/*) = Trs ( X] ® ^ ® ® Ri)*) • 

2 = 1 


By linearity, the previous equality gives 


VT GAfnfc(C), TrB(17-y-C/*) = TrB( 

2=1 


Then, by definition of the partial trace, we obtain the following equivalences 

r 

vy G Mnk(C), A G Mn(C), Tr(U ■ Y ■ U*A ® I) = Tr(^ ^ ® i?* • y • ([/* ® Ri)*A ® l) 

2 = 1 
r 

vy G Mnk(C), A G Mn(C}, Tr(Y ■ U*A (g) lU) = Tt (y ■ U*AUi ® Mi) 

2=1 

r 

VA^MniC), U*A®IU = YRiXUi®Mi. ( 35 ) 

2 = 1 


Let us now study (35). Consider two matrices A, B in Mn{C). Since U*AB (g) lU = U*A (g) lU ■ 
U*B ® lU, we obtain from (35) 

r r 

Y U-ABUi ®Mi=Y U*AUiU*BUj ® MiMj. (36) 

i=l *J=1 

Applying (36) with A = ee* and B = ff* where e and / are orthogonal vectors of C", it directly 
follows 

0 = ^ U*ee*UiU*ff*Uj ® MiMj. 
i¥=j 

Taking the trace, we obtain 

Ve T /, 0 = Y me,U;f)\^ Tt{M,M,) = Y \{e,UiU*f)\^ Tr(MiM,-). 

i¥=j i¥=j 

Note that the previous equation is actually a sum of non-negative terms equals to 0. Therefore, we 
conclude that for i ^ j 


TT{MiM,) = 0 or VeT/, {e,UiU*f) = 0. (37) 

Now since for all i ^ j, Ui Uj, the unitary matrix UiU* is not a multiple of the identity. Thus, 
we claim that we can find orthogonal vectors e and /, such that (e, UiUj f) ^ 0. Indeed, the matrix 
UiU* is diagonalizable in an orthonormal basis (rtp)p=i,..,,n with related eigenvalues {pp)p=i,...^n- 
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Since UiU* 0 Cl, all the eigenvalues are not equal, for instance / H 2 - Let us consider the 
orthogonal vectors e = U 1 — U 2 and f = Ui+U 2 - We can easily check that (e, UiU* f) = —/r 2 / 0. 

Finally, we deduce from (37) that Tr(MjMj) = 0 for all i ^ j and thus MiMj = 0. The Mj’s 
being orthogonal and sum to the identity, we define the matrix V = Yli=i Li < 8 ) Ri belonging to U^. 
We can now easily check that the unitary operators U and V induce the same quantum channels 
for all /3 in By Lemma 1.3, there exists W of Uk such that U = {I ® W)V] this proves 

the result. □ 


Proposition 7.2. If n = 2, then ^^ocfc-diag ~ Uunitai- In particular, the chain of inclusions (34) 
collapses: 


U, 


block—diag Upj,Q},i—iin Upj-Qj, U^yiixed U^inital- 


Proof. Consider { 61 , 62 } the canonical orthonormal basis of and consider a matrix U in Uunitai 
written as 2 x 2 block matrices 


U = 


A B 
C D 


eie\ (8) A + 6162 ® B + 626} (8 C + 6263 ® D 


with A, B, C and D in Aik{C). As proved in Theorem 3.1, both U and (U^)'^ are unitary matrices 
and therefore 

AA* + BB* = A* A + B*B = I, AA* + CC* = A* A + C*C = I, 

CC* + DD* = C*C + D*D = I, DD* + BB* = D*D + B*B = I, (38) 

AC* + BD* = A*C + B*D = 0, AB* + CD* = A*B + C*D = 0. 


Our aim is to applied Proposition 6.4. Consider the orthonormal basis {Eij = eje*,i,j = 1,2} of 
A 42 (C), it is clear that 

{XijX*,i,i,j, k,l = 1,2} = {XY* \X,Y = A, B, C, D} 

{X*jXki,i,j, k,l = 1,2} = (XW \X,Y = A, B, C, D}. 

Let us now prove that {XY* \ X,Y = A,B, C, D} and {X*Y \X,Y = A, B, C, D} consist of normal 
commuting matrices. It can be noticed that it is sufficient to check that, for all X,Y,Z = A, B,C, D, 

XY*Z = ZY*X. 


Using the symmetry A, D in (38) together with the symmetry B, C, the 64 different cases boil down 
to 11 non-trivial cases : AA*B,AA*D,AB*B,AB*C,AB*D,AD*B,AD*D,BA*C,BA*D,BB*C 
and BC*D. Each of them can be easily checked as for instance 

AA*B = {I- BB*)B = B{I - B*B) = BA*A 


or 


BC*D = B{C*D) = B{-A*B) = {-BA*)B = DC*B. 

The result then holds. □ 

Remark 7.3. A similar result has been obtained in [14, Theorem 9], under more stringent as¬ 
sumptions. More precisely, it is shown in [14] that, when n = 2, ~ Uprob, assuming 

that the unitary operators appearing in the mixed-unitary decomposition of channels are linearly 
independent. 

Swapping the roles of n and k, we obtain the following result. 

Proposition 7.4. If k = 2, then U^ock-diag = I^unitai ■ 


ON SOME CLASSES OF BIPARTITE UNITARY OPERATORS 


21 


8. Conclusions and open questions 

We end this work with a list of questions that we have left unanswered (or even untouched). We 
hope to get back to some of these problems in some future work. 

We start with the problem of computing the dimension of the algebraic variety Uunitai] recall 
that previously, we have looked at the enveloping tangent space of this variety, at some particular 
points. 

Conjecture 8.1. Show that dimUunitai = kn‘^ + nk‘^ — nk. 

It has been showed in Theorem 7.1 that any operator in the set lAprob-Un (which is a subset of 
Umixed) is block diagonal, with respect to the system A. Moreover, in the qubit case n = 2, we 
have ~ ^mixedi See Proposition 7.2. We conjecture that this equality always hold, and 

that the technical restrictions appearing in the definition of Uprob-Un are actually superfluous. 

Conjecture 8.2. For all values ofn,k, it holds that 

^block-diag ~ ^mixed- 

Regarding bipartite unitary operators producing PPT channels, we have left the following prob¬ 
lem open. 

Question 8.3. Given a subspace R C C" (8) C^, let P e Al^'^(C) be the orthogonal projection on 
V. Characterize the set of subspaces V such that 

Al*1fc(C) 9 (4 ® P){Fn ® 4)(4 ® R) > 0, 

where F is the flip operator. 

This brings us to the problem of characterizing the set Ueb and comparing it to Uppt (at the 
level of quantum states, this would be the fact that the PPT criterion for separability is necessary 
in all dimensions, and sufficient for nk < 6). 

Question 8.4. Provide a description of the set IAeb- For which values of n,k, is it true that 
UppT = Ueb ? 

At the level of examples, beside the obvious inclusion Uconst P Geb, we also have^, when n = k, 

^block—diag ’ — ^EB- 

Indeed, for a unitary operator U = (X]r=i ® ^ifi) '^n, the corresponding quantum channel reads 

n 

i=l 

which is entanglement-breaking. 

Finally, we consider the following classe of bipartite unitary matrices yielding channels of interest 
in quantum information theory. 

Gcq = {U G ^7 I V/3 G A4],’''~(C), is a classical-quantum channel} 

Gqc = {U G 77 I V/3 G A7],’''~(C), Lu,j 3 is a quantum-classical channel} 

77cc = {U G 77 I V/3 G A7].’''~(C), Lp^p is a classical-classical channel} 

The study of these classes has been initiated in [14, 15], where mainly the qubit case n = 2 has 
been discussed. The structure of these operators in the general case remains open. 

Question 8.5. Characterize the setsUcQ, Uqc, and Idee- 

^We thank Siddharth Karumanchi for pointing this out to us. 
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Appendix A. Necessary and sufficient conditions for the existence of a block-SVD 

In this section, we establish necessary and sufficient conditions for the existence of a block singular 
value decomposition of a bipartite operator with respect to the second sub-system B. Moreover, 
we present an algorithm for obtaining such a decomposition when it does exist. These results are 
inspired from [3, 7], see also [12, Theorem (2.5.5) and Section 7.3, Problem 25]. 

We shall denote, for any matrices X G Wf„(C) (8) AikiC) and A G Wf„(C), 

Aa := [T7®id](A-A® 4 ). 

Theorem A.l. Let X G Af„(C) (8 Alfc(C). The following conditions are equivalent: 

(1) The matrix X has a block-diagonal SVD; there exists matrices Xi,...Xp G Afn(C) and 
partial isometries Ri,...,Rp G Affc(C) having orthogonal initial, resp. final, projections 
such that 

p 

X = '^Xi® Ri. 

i=l 

(2) The families f^sp. {Xj^XB}A,BeMn(C)> consist of commuting, normal 

operators. 

2 2 

(3) The families resp. consist of commuting, normal operators. 

We write Xa ■= Xe* for an orthonormal basis {Ea}a=i of M-n{C), in such a way that 

X = Eii E^®x^. 

We denote the set of matrices X satisfying the above condition(s) by Ai^ock-diag 

k 

Mtck-^ag ■■= {X e MnkiC) | A = A, ® ej*}. 

i=l 

The set A^^ocfc-diag defined in a similar way, by swapping the roles of the two tensor factors. 

Proof. The implication (2) (3) is obvious. Let us first show (1) (2). For a matrix A as 

in (32), we have 

p 

XaX% = J^Tr(A,A)Tr(A,B)Pi, 

i=l 

where Pi = RiR* is the orthogonal projection on the image of the partial isometry Ri. This relation 
shows that the matrices {XAXfijA,B&Mn{C) normal and diagonalizable in the same orthonormal 
basis of C". A similar argument shows the result for the matrices {X\Xb}A,BGM n{C) 

Let us now show (3) (1). The fact that the normal matrices AqA^ commute implies they 

have the same set of eigenprojectors P*: 

p 

Va,/3, A,A^ = ^AS“’^)Pi. (39) 

i=l 

In the same vein, we have, for another set of orthogonal eigenprojectors Qf. 

p' 

XlXp = Y,Tt’^^Qi- ( 40 ) 

i=\ 

Letting a = /? in (39) and (40) and using the fact that the matrices A^ A* and A* A^ have the same 
(positive) eigenvalues (counting multiplicities), we have that p = p' and there exists permutations 
ctq, G Sp and complex numbers such that 
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for some partial isometries having initial projection Qaa{i) final projection Pi. Plugging 
the last expression into (39) and (40), we find that the permutations cTq, must be equal; we shall 
assume, by re-ordering the eigenprojectors Qi, that these permutations are all equal to the identity. 
Using similar arguments, the partial isometries cannot depend on a, and we write R-'^^ = Ri- 
We have thus 

i=l 


and thus 


Setting 


concludes the proof. 


X = J2Ea^ ^ 

cx=l \i=l 

E ( EE’-E” I ®-«‘' 

i=\ \a=l 


Y := E 


Q. = l 


□ 


Remark A. 2. The above result provides us with an efficient way of checking whether a given 
bipartite operator X has a block-SVD: pick a basis of XiniC) (e.g. the usual matrix units) and 
check the condition in (3) above. 


Remark A. 3. Restricting to a = 13 in the condition (3) does not yield an equivalent statement, as 
it is shown in the example below. For X G A4 2 . 2 (C) given by 


X = 


Ai 2 — A21 — 


(41) 


with the choice of the canonical matrix units for A42(C), the matrices Xa are Xu = 


, and X 22 = l 2 - All the matrices X^X^ are diagonal, hence they are normal 


and commute. However, the matrices X 11 X 22 and X 12 X 22 do not commute, so the 4x4 matrix X 
does not satisfy the equivalent conditions from Theorem A.l, hence it does not have a block-SVD. 


Corollary A. 4. The set is a real algebraic variety. 


Proof. A matrix X belongs to if and only if the two families of matrices {AqA^} 

and {A*A^} commute; these commutations conditions can be restated as (degree 4) polynomial 
conditions in the real and imaginary parts of the elements of A. □ 


Before computing in the next proposition the dimension of the real algebraic variety 
we would like to give a heuristic argument in the form of parameter counting. The choice of the 
two orthonormal bases {cj} and {/*} in (32) corresponds to a total of 2k‘^ real parameters, since 
diiUR^fc = Each matrix Aj accounts for 2n^ real parameters, so, in total, we get 2kn^ extra 
real parameters. However, in each term Aj ® eif*, two of the three complex phases of Aj, e^, /* are 
redundant, so we have over counted 2k‘^ real parameters. We conclude that the real dimension of 
■A^Uock-diag should be 2k‘^ + 2kn‘^ — 2k, fact which we rigorously prove next. 
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Proposition A. 5. The real dimension of the algebraic variety -M-biock-diag 2fc(n^ + A: — 1). 

Proof. For the terminology and the results used in this proof, we refer the reader to [11, Chapter 
11]. Let us introduce the flag manifold (see [11, Example 8.34] or [5, Section 4.9]) 

Fh=Uk/{U^), 

which has real dimension dim® Fl^ = — k. Consider the map 

V; : Mn{Cf X Fife X FU ^ M^ck —diag ^ Mn(C) 0 Mk(C) 

k 

E A ® ei/;- 

i=l 

Obviously, ip is surjective, so we have < 2k{n^ + A — 1). To get the reverse 

inequality, define A4n(C)^ to be the set of pairwise distinct A-tuples of matrices. It is trivial 

to check that the map p, obtained by restricting p to J\An(C)^, is /c!-to-one, so the conclusion 
follows. □ 
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